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1. Introduction 

Testing the validity of Chiral Perturbation Theory (ChPT) is a challenging task because of the 
many unknown parameters, called the Low Energy Constants (LECs), entering into the theory. In 
particular, at NNLO 90 unknown constants, the Q, appear in the Lagrangian. 

One way to overcome this problem is to study different combinations of observables that 
depend on the C, in the same way. These lead to C, -independent relations which can be used to 
perform the test. Furthermore those combinations might be useful to gain information on the LECs 
too, since they let us isolate the same combinations of C, using different observables. 

In we study 76 observables at NNLO and find 36 such relations. We compare ChPT NNLO 
predictions with data/dispersive results for 13 of these. The observables involved are the ones in nn 
and ;r^r-scattering and in Ki4 decay. Here we first discuss how we perform the numerical analysis, 
the results of which appear in Tab.|I|, ^ |3| and ^ then we present for each process the relations 
studied. Finally we show some preliminary results for a new global fit of the L, at NNLO. 

2. Numerical Analysis 

The numerical analysis of the C, -independent relations has been done in the following way. 
First we evaluate the combinations of observables appearing in each side of the relations using 
experiment/dispersive (exp) results of for nn scattering, for kK scattering and ^ |6|] for 
K(4 decay. Then we use ChPT results up to order [^, |8|, ^] setting the L, to the values of fit 10 
in [|l|]. Finally we subtract from the first (exp) evaluation the ChPT one. These differences will 
contain the C, part and higher order corrections. They have been quoted in Tab.|l|, ^, ^ and ^ in 
the columns labeled remainder. This has been done for each side of the relations under study. To 
check whether a relation is well satisfied we compare the remainders of its left-hand-side (LHS) 
and right-hand-side (RHS). Since they contain the same Q combinations, they should be equal 
within the uncertainties. 

The errors quoted in the second columns of Tab.|], ^, ^ and ^ aie obtained adding in quadrature 
the uncertainties in [|3[ ^, ^ This might result in an underestimate of the total error because 
of correlations. The theoretical errors due to the NLO LECs are shown in brackets in the columns 
of Tab.|I], ^, ^ and ^ labeled NNLO 1-loop. They are obtained by varying all the L, around the 
central values of fit 10 according to the full covariance matrix as obtained by the authors of [|l]] 
and exploring the region with ;t;^/dof « 1 . The error is then estimated as the maximum deviation 
observed. The error for the L, contribution at NLO is never shown since it drops out of all the 
relations. No uncertainties due to higher order contributions have been added. The uncertainties 
due to theoretical errors are mostly on the last quoted digit. 

3. TtTt scattering 

The Tin scattering amplitude can be written as a function A{s,t,u) which is symmetric in t,u: 
A{Tt"K^ ^ k'k'^) = 5"'''5''''A{s,t,u) + 5''''d'''''A{t,u,s) + d"''' d''''A{u,t,s) , (3.1) 
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where s,t,u are the usual Mandelstam variables. The isospin amplitudes T\s,t) (/ = 0, 1,2) are 
T^{s,t) = 3A(s,f,M) +A{t,u,s) +A{u,s,t), T\s,t) =A{s,t,u) -A{u,s,t), T^{s,t) =A{t,u,s) + 
A{u,s,t) , and are expanded in partial waves 

T\s,t) = 327r£(2^+l)P/;(cose)4(5), (3.2) 

1=0 

where t and u have been written d&t = —^{s — Am\){\ —cos d), u = — 2(5 — 4m|)(l +cos d). Near 
threshold the t\ are further expanded in terms of the threshold parameters 

tl{s)=q^\a[+b\q^ + i?{q'')), = \{s-Aml), (3.3) 

where a\,b\ ... are the scattering lengths, slopes, We studied the 11 parameters where a depen- 
dence on the Cj shows up. Using s + t + u = 4m\ we can write the amplitude to order as 

A{s^t^u) = b\ + b2S + b2,s^ + b^{t — uf' + b=,s^ + b(,s{t — m)^ + non polynomial part (3.4) 

The tree level Feynman diagrams give polynomial contributions to A{s,t,u) which must be ex- 
pressible in terms of bi, . . . ,b(,. Therefore we expect and find 5 relations: 

[5Z7g-2Z?[]-27a}-15ag + 6a{j]^. = -n[b\]^_, (3.5) 

[3a\ + bl] ^_ = 20 [bj -b^2-4 + 41 c, , (3-6) 

[b^o + 5bl + 9a\]^_ = 90 [4 - b'i] ^_ , (3.7) 

[3Z7} + 25a^]^_ = 10[a^]^., (3.8) 

[Sbj + lb'^]^^ = 21 [al]^^, (3.9) 

where [A]q = C[-dependent part of A. All quantities are expressed in units of ni^+- In fact, since 
these relations hold for every contribution to the polynomial part, they are valid for the NLO tree 
level contribution as well and for two- and three-flavour ChPT. Thus they get L, -contributions only 



at NNLO via the non polynomial part of Eq. (3.4) 



In Tab. [T] we show our numerical results. We quote the left-hand-side (LHS) and right-hand- 
side (RHS) of each of the relations. In the second column we use the values of the threshold 
parameters of The next columns use the ChPT results of ^ and give the contributions from 
pure one-loop at NLO, the tree level NLO contribution, the pure two-loop contribution, and the 
dependent part at NNLO (called NNLO 1-loop). 

Comparing the remainders of the LHS with the RHS ones, we see that the first three relations 
are very well satisfied, while the last two work at a level around two sigma. 

We can also check how the two-flavour predictions hold up. Since here the corrections are in 
powers of rather than in powers of m^, the expansion should converge better. For the ChPT 
evaluation we use the threshold parameters as quoted in ^ for their best fit of the NLO LECs. 
The result is shown in Tab. |^ We see the same pattern as for the three-flavour case: the first three 
relations are very well satisfied while the last two are somewhat worse but below two sigma. 
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LHS (O) 



RHS (3.5) 



10 LHS (3.6) 



10 RHS (3.6) 



LHS (^ 
RHS (0) 
10 LHS (|3j) 
10 RHS (U) 
10^ LHS (0) 
10^ RHS (^ 



i 



0.009 ±0.039 
-0.102 ±0.002 
0.334 ±0.019 
0.322 ±0.008 
0.216±0.010 
0.189 ±0.003 
0.213 ±0.005 
0.175 ±0.003 
0.92 ±0.07 
1.18±0.04 



NLO 


NLO 


NNLO 


NNLO 


remainder 


1-loop 


LLCs 


2-loop 


1-loop 




0.054 


-0.044 


-0.041 


-0.002(3) 


0.041 ±0.039 


-0.009 


-0.044 


-0.060 


-0.008(6) 


0.018 ±0.002 


0.209 


0.097 


0.103 


0.029(11) 


-0.105 ±0.019 


0.177 


0.097 


0.120 


0.034(13) 


-0.107 ±0.008 


0.166 


0.029 


0.053 


0.016(6) 


-0.047 ±0.010 


0.145 


0.029 


0.049 


0.020(7) 


-0.054 ±0.003 


0.137 


0.032 


0.053 


0.035(12) 


-0.043 ±0.005 


0.121 


0.032 


0.050 


0.029(10) 


-0.057 ±0.003 


0.36 


0.00 


0.56 


-0.01(13) 


0.00 ±0.07 


0.42 


0.00 


0.57 


0.03(13) 


0.15 ±0.04 



Table 1: The relations found in the TTTT-scattering. The lowest order contribution is always zero by con- 
struction. The NLO LEC part satisfies the relation, as it should. Notice the extra factors of ten for some of 
them. All quantities are in the units of powers of m.,^+ . 



LHS (3.5) 



RHS (3.5) 



10 LHS (3.6) 



10 RHS (3.6) 



LHS (^) 

RHS 
10 LHS ^ 
10 RHS dJ) 
10^ LHS (U) 
10^ RHS (^ 



i 



0.009 ±0.039 
-0.102 ±0.002 
0.334 ±0.019 
0.322 ±0.008 
0.216±0.010 
0.189 ±0.003 
0.213 ±0.005 
0.175 ±0.003 
0.92 ±0.07 
1.18±0.04 



two-flavour 



-0.003 
-0.097 
0.332 
0.318 
0.206 
0.189 
0.204 
0.176 
1.00 
1.15 



remainder 



0.007 ±0.039 
-0.005 ±0.002 
0.002 ±0.019 
0.004 ±0.075 
0.010 ±0.010 
0.000 ±0.003 
0.009 ±0.005 
-0.001 ±0.003 
-0.08 ±0.07 
0.04 ±0.04 



Table 2: The relations found in the TTTT-scattering evaluated in two-flavour ChPT. In the second column we 
have used the NNLO results quoted in [^. Notice the extra factors of ten for some of them. All quantities 
are in units of powers of m;^+ . 



4. :;r^ scattering 

The nK scattering has amplitudes T\s,t,u) in the isospin channels / = 1/2,3/2. As for Tin 
scattering we introduce the partial wave expansion of the isospin amplitudes 



-|-oo 



and we define scattering lengths a^, Z?' by expanding the f'(i') near threshold: 



(4.1) 
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and t = —lq\fr{ \ — cos 0), u = —s — t + 2m\ + 2m|. Again we studied only those observables 
where a dependence on the Q shows up. 

It is also customary to introduce the crossing symmetric and antisymmetric amplitudes {s,t, u) 

3T+{s,t,u) = T^l'^{s,t,u) + T^l^{s,t,u), {s,t,u) = T^^^{s,t,u) -T^/^{s,t,u), (4.2) 

wliich can be expanded around t = 0, s = u using v = (s — u)/ (4m^) (subtlireshold expansion): 

OO CO 

T+is,t,u) = X cp'v^J, T-{s,t,u) = X cTjt'v^^+K (4.3) 
ij=0 i,j=0 

There are 10 subthreshold parameters that have tree level contributions from the NNLO LECs. In 
Cqj and the same combination — Ci +2C3 + 2C4 appears [||]: 

16p2[c2o]c, = 3[coi]c,. (4.4) 

Therefore in the isospin odd channel only three subthreshold parameters get independent con- 
tributions from the Q. So for the 7 differences aj = — a^/^ and bj = b^^"^ — Z?^^^ getting 
contributions at NNLO and three subthreshold parameters we expect four relations: 



(p^ + 3p3 + 3p + 1) [ar] c, = 2p' (P + 1)' {K\ ^_ - ^p (p2 + 1) c, 

+ + + (4.5) 

5 (P + 1)^ hie, = K]c, - ^^^-^ K1c, + ^^^i^ [^0 ]c, > (4-6) 

5 (p2 + 1) [a^ ]c, = [«r]c, + 2P [^r]c, ' (4.7) 
7 (P' + 1) [«3lc, = h1c, +2P [^2 ]c,- ' (4-8) 

the threshold parameters are expressed in units of m-,^^ and we use the symbol p = niK/mj^. 

1 /2 3 /2 1 /2 

r+ brings in 7 more combinations of threshold parameters, a'^ =a/ + 2a ^ and bf =b/ + 

3 /2 

2^^' , but there are 6 independent subthreshold parameters so we find only one more relation: 

7 [4] c, = Tp ["^] c. - H] c, + ] c, - 6^ [4] c. - 3^ [^o1 c, • (4-9) 

Again these relations hold for all tree-level contributions up to NNLO. The numerical check 
is shown in Tab. ||. The columns in Tab. ^ have the same meaning as in Tab. 

The first relation is reasonably satisfied, somewhat below two sigma. The second relation has 
a large discrepancy but if we assume a theory error of about half the NNLO contribution it seems 
reasonable. The third relation is well satisfied but the RHS has a rather large experimental error. 
The fourth relation does not work well, mainly due to the fact that we seem to underestimate the 
value for . The last relation works well. 
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NLO 
1-loop 



NLO 
LLCs 



NNLO 
2-loop 



NNLO 
1-loop 



remainder 



LHS Q 
RHSQ 
10 LHS Q 
10 RHS Q 



100 LHS (4.6) 



100 RHS (4.6) 



100 LHS (4.8) 



100 RHS (^) 

10^ LHS 

10^ RHS (|49|) 



5.4 ±0.3 
6.9 ±0.6 
0.32 ±0.01 
0.37 ±0.01 
-0.49 ±0.02 
-0.85 ±0.60 
0.13±0.01 
0.01 ±0.01 
0.29 ±0.03 
0.31 ±0.07 



0.16 
0.42 
0.03 
0.02 
0.08 
0.03 
0.04 
0.01 
0.09 
0.03 



0.97 
0.97 
0.12 
0.12 
-0.25 
-0.25 
0.00 
0.00 
0.00 
0.00 



0.77 
0.77 
0.11 
0.10 
-0.17 
0.11 
0.01 
0.00 
0.06 
0.06 



-0.11(11 

-0.03(7 
0.00(2 

-0.01(2 
0.05(3 

-0.03(13 
0.03(1 
0.00(1 
0.01(2 
0.05(3 



0.6 ±0.3 
1.8±0.6 
0.07 ±0.01 
0.14±0.01 
-0.21 ±0.02 
-0.71 ±0.60 
0.05 ±0.01 
-0.01 ±0.01 
0.13 ±0.03 
0.17 ±0.07 



Table 3: The relations found in the TT/T-scattering. The tree level contribution to the LHS and RHS of 
relation 1 is 3.01 and vanishes for the others. The NLO LECs part satisfies the relation. Notice the extra 
factors of ten for some of them. All quantities are in the units of powers of m„+ 

5. Ttn and nK scattering 

If we consider the nn and nK system together we get two more relations due to the identities 



[^5]c, = [4o]c+-[' 



20j C, 



16p2 



(5.1) 



where ■ {c'^-) are expressed in units of m^'^^^^\m^'^^^). We can express these relations in terms 



of the threshold parameters (all quantities expressed in powers of 171,^+): 



:i+p) 



7(l-4p±p2) [aj]^^+[a+ + 2pb. 



2 Jc, 



(5.2) 
. (5.3) 



The numerical results are quoted in Tab. ^. The first relation does not work but the second is well 
satisfied. If we look in the numerical results we see that plays a minor role in the RHS of the 



second relation but is important in the first, so this could be the same problem of relation (4.8). A 
related analysis can be found in [p^. 



6. K, 



£4 



The decay K^{p) n^{p])n {p2)e^{pe)v{pv) is given by the amplitude [|TTj] 



^v:APvW-r5)viPi)iv^-A^) 



(6.1) 



where and are parametrized in terms of four formfactors: F, G, H and R (but the R- 
formfactor is negligible in decays with an electron in the final state). Using partial wave expansion 



and neglecting d wave terms one obtains [ |I2| ] : 

F = /. ±/,V ±/:V ±/>./4m2 ±/,(7;,Xcos d + ... 
Gp = gp+ g'q^ + g'q^ + g'^Se/Am\+ gtO^Xco^ B + ... 



(6.2) 
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NLO 
1-loop 


NLO 
LLCs 


NNLO 
2-loop 


NNLO 
Lloop 


remainder 


10^ LHS ( 
10^ RHS ( 
10 LHS ( 
10 RHS ( 


5.2 
5.2 

5.3|; 

5.3 


) 
) 

) 


0.34 ±0.01 
0.38 ±0.03 
-0.13±0.01 
-0.09 ±0.02 


0.12 
0.12 
-0.12 
-0.05 


0.00 
0.00 
0.00 
0.00 


0.16 
0.05 
-0.05 
-0.02 


0.00(4) 
0.04(2) 
0.02(2) 
-0.01(1) 


0.05 ±0.01 
0.16±0.03 
0.01 ±0.01 
-0.01 ±0.02 


LHS (6.4) 
RHS (6.4) 


-0.73 ±0.10 
0.50 ±0.07 


-0.23 
0.19 


0.00 
0.00 


-0.15 
0.10 


-0.05(6) 
0.03(4) 


-0.29 ±0.10 
0.18 ±0.07 



Table 4: The relations found between nn and Tr/T-scattering lengths and between the curvature in F in /r^4 
and nK scattering. All quantities are in the units of powers of mj^+ . 



Here s^{se) is the invariant mass of dipion (dilepton) system, and 
of the pion in their rest frame w.r.t. the kaon momentum and t — u 
ChPT results [|8|, ^ we find one relation between the quantities defined in (5J) and nK scattering 



Sj^/ (4m|) — 1. 6 is the angle 
-iGj^Xcosd. Using NNLO 



V2[f':]^=64pF4c 



30j C, 



(6.3) 



This leads to a relation between nK threshold parameters and which, with all quantities ex- 
pressed in units of m;j+, reads: 



v2[/r: 



c, 



3271- 



1±P 



^(2 + P±2p^) 



2 Jq 



(6.4) 



Numerical results for ( 5.4) are shown in Tab. ^. The experimental results is taken from |j5|] for 
fs /fs and from for /j. This should be an acceptable combination since the central value for 
f'Jfs and //s from are in good agreement with those of [Q]. This relation is not satisfied: the 
sign is even different on the two sides. Notice that, in both cases, we also see that the ChPT series 
has a large NNLO contribution. 

It has been already noticed, see [|]] and Fig. |l|, that ChPT, at present, underestimates the cur- 
vature f". On the other hand there are indications that dispersive analysis techniques might help 
solving this problem: Fig. 7 in [[T]] shows that the dispersive result of [13] has a larger curvature 
then the two-loop result. Therefore, we do not consider this discrepancy a major problem for ChPT. 



7. New fits of the NLO constants (preliminary results) 

As remarked in [[l4]], many NNLO calculations are now available in three-flavour ChPT. Be- 
sides, new lattice and dispersive results and further experimental data are at our disposal too. A 
study of the predictive power of NNLO ChPT is needed, and therefore also an update of the L, fit. 
For this reason we are working on a new program to perform this fit with many more observables 
implemented. So far we have included masses and decay constants, K£4 formfactors, nn and nK 
scattering lengths and the scalar pion radius. For now we rely on the resonance estimates of the Q 
used in although our plan is to achieve more information on them. 

Our first preliminary results are summarized in Tab. ^ In the second column we quote fit 10 
of This was found using the available linear fit for K£4 of [^, Fk/Fji = 1.22, the kaon and 
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6 

5 



4 



CO n 



2 



1 





0.02 0.04 0.06 0.08 0.1 0.12 0.14 0.16 

7 

Figure 1: The absolute value of the formfactor at S£ ~ cos = as a function of (in Gev units) 
above and below threshold. The NNLO result nicely reproduces the linear fit quoted in 1^, but not the large 
negative curvature in [|[]. The line at the bottom is the contribution coming from the C,-, which has a positive 
curvature. 



fit 10 [jl]l fit 10 iso NA48/2 Fk/F^ All 





0.43 


0.40±0.12 




0.98 




0.97 


0.99 ±0.13 


lO^L^ 


0.73 


0.76±0.12 




0.78 




0.79 


0.60 ±0.22 


lO^L^ 


-2.35 


-2.40 ±0.37 


-3.14 


-3.12 


-3.07 ±0.59 


lO^L^ 


= 


= 




= 




= 


0.65 ±0.64 


lO^L^ 


0.97 


0.97 ±0.11 




0.93 




0.72 




0.53±0.10 




= 


= 




= 




= 


0.07 ±0.65 


lO^L^ 


-0.31 


-0.30 ±0.15 




-0.30 


-0.26 


-0.21 ±0.15 


lO^L^ 


0.59 


0.61 ±0.20 




0.59 




0.48 


0.37 ±0.17 


(dof) 




0.25 (1) 


0.17 (1) 


0.19 (1) 


0.78 (4) 



Table 5: Preliminary results for the fits. Lg = 0.59 x 10 ^ everywhere, as found from the vector pion radius 
in (jTsjl. See text for a longer discussion 




LO 

NLO 
NNLO 
reso only 
NA48 
E865 linear 



eta masses with isospin breaking corrections included and setting L4 = Le = 0. In the column 
labeled fit 10 iso we quote the fit we find using the same input as fit 10 but without including 
isospin breaking. As you see the two fits are in good agreement. The column NA48/2 relies on the 
new experimental data from [^]. We checked that the fit does not change including the curvature 
With this fit ChPT predicts the value f'J = —0.90 to be compared with the experimental one 
f'J = —1.58 ±0.064. Note that the fit in ^ shows large correlations between the slope and the 
curvature of the F^ formfactor which have not been taken into account yet. The values of Li and L3 
change drastically. The third column shows the fit obtained changing the ratio Fk/F-^ to 1.19. This 
affects mainly L5 and L^. The last column shows the fit obtained letting L4 and Lg free, and adding 
ajj, cxq, aj^ , Oq^^ and the scalar pion radius. The value obtained for L4 is larger then expected. 



Some more comment can be found in [ 16] 
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8. Conclusions 

We have performed a systematic search for relations between observables that allow a test of 
ChPT at NNLO order in a Q-independent way. We studied in detail the relations for the nn, nK 
scattering and Ke^ since for these cases enough experimental and/or dispersion theory results exist. 

The resulting picture is that ChPT at NNLO mostly works but there are troublesome cases. 
The KK system alone works well. The nK system alone works satisfactorily but with some dis- 
crepancies. The same can be said for the combinations of both systems. A common part in these 
two cases is the presence of . Comparing nK scattering and Ki4 leads to a clear contradiction 
which needs further investigation. 
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